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ABSTRACT 



We present an empirical near-infrared (NIR) spectroscopic method for estimating M 
dwarf metalhcities, based on features in the ff-band, as well as an implementation of a 
similar published method in the K-haiid. We obtained i?~2000 NIR spectra of a sample 
of M dwarfs using the NASA IRTF-SpeX spectrograph, including 22 M dwarf metallicity 
calibration targets that have FGK companions with known metalhcities. The H-hand 
and K-band calibrations provide equivalent fits to the metalhcities of these binaries, 
with an accuracy of ±0.12 dex. We derive the first empirically calibrated spectroscopic 
metallicity estimate for the giant planet-hosting M dwarf GJ 317, confirming its super- 
solar metallicity. Combining this result with observations of eight other M dwarf planet 
hosts, we find that M dwarfs with giant planets are preferentially metal-rich compared 
to those that host less massive planets. Our i7-band calibration relies on strongly 
metallicity- dependent features in the i?-band, which will be useful in compositional 
studies using mid to high resolution NIR M dwarf spectra, such as those produced by 
multiplexed surveys like SDSS-III APOGEE. These results will also be immediately 
useful for ongoing spectroscopic surveys of M dwarfs. 



Subject headings: planetary systems 
niques: spectroscopic 



stars: abundances — stars: late-type — tech- 



Introduction 



Chemical analysis of M dwarfs is notoriously difficu lt due to the complex molecular spectra 



of their cool atmospheres (jMould 



1976 



Gustafssonlll989l ). Substantial progress in metallicity es- 



timation for M dwarfs has been made recently, both through photometry and spectroscopy, and 
motivated i n large part by by the question of whether the metallicity- giant planet (> 0.5Mjupitcr) 



correlation ( Gonzalez 



1997 



Santos et al 



2004 



Fischer &: Valentill2005l ) extends to low-mass stars. 
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Photometric techniques are b ased on the higher o pt ical opacity in meta l-rich stars, which shifts 
visible light into the infrared (e.g. iDelfosse et al.ll200Cll ). iBonfils et al.l (|2005l ) first calibrated this re- 
lationship, and later groups h av e modified it and a ddressed its limitations (jjohnson fc Appsll2009l : 



Schlaufman Sz Lau ghliJ 2010). iNeves et al.l (l201lh evaluated the published photometric calibra- 
tions, and found Sch laufman &: Laughlinl (|2010l ) to have the lowest dispersion (rms = 0.19 it 0.03 



dex) . Photometric techniques have generally bee n limited by requiring absolute magnitudes of the 
targets, although the new technique presented in iJohnson et al.l (120111 ) depends on color alone. 



Empiri cal spectroscopic technique s have also been developed, whic h rely on calibrations us- 



ing optical (jWoolf &: WallersteinI |2006| ) or near-infrared (NIR) ii'-band (jRoias-Avala et al.l l20ld . 
hereafter, RAlO) feature strengths. RAIO used Na i, Ca i, and H2O features to model metallic- 
ity, and calibrated their model using M dwarf companions to FGK stars, with an uncertainty of 



(Bean et al. 


2006: 


Maness et al. 


2007; 


Onehag et al. 


2011 



difficult to be used on a large scale due to the intrinsic faintness of M dwarfs, their dense spectra, 
and incomplete linelists. 



We are developing a target list in support of the Habitable Zone Planet Finder (iMahadevan et al. 



2OIOI ). a NIR precision radial velocity spectrograph for finding planets around M dwarfs. Metallicity 
estimates will be an important parameter of this target list, so that we may sample a wide range 
of stellar metallicities. Ou r overall list is drawn from the all-sky nearby M dwarf catalog presented 



m 



Lepine Sz GaidosI (1201 ll ). and we are carrying out a campaign on the NASA Infrared Telescope 
Facility (IRTF) to obtain i?~2000 spectra of these stars and estimate their metallicities using the 
method described herein. M dwarfs will also be ta rgeted in high-res olution surveys, such as th e 



Sloan Digital Sky Survey III (jEisenstein et al 



2011 



APOGEE survey (lAllende Prieto et al.ll2008l ') 



which will deliver i2~20,000 H-hand spectra of hundreds of M dwarfs at multiple epochs. 

Here, we report on a new NIR spectroscopic calibration of metallicity-sensitive features in 
the -ff-band, calibrated with well-separated M dwarf companions to FGK stars with known metal- 
licities. This calibration yields a smaller scatter than the RAIO K-hand calibration. The main 
promise of the /f-band technique is that the relevant metallicity-sensitive features will be included 
in APOGEE high-resolution H-hand spectra of M dwarfs, allowing a high-resolution metallicity 
calibration and potentially direct spectroscopic modeling like that accomplished by Oil. In §2, 
we describe our observations, in §3, we explain our analysis and error estimates, and in §4, we 
discuss our calibrations and our analysis of M dwarf planet host metallicities, which indicates that 
M dwarfs with giant planets are preferentially metal-rich compared to those that host less massive 
planets. 



■^At submission of this paper, 
with an uncertainty of 0.14 dex. 
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Observations 



During A ugust and Novemb er 2011, we observed more than 200 M dwarfs using the SpeX 
spectrograph (jRavner et al.ll2003l ) on the NASA IRTF. For cahbration of our empirical metahicity 
scale, our target list included 22 M dwarfs in binary sy stems with FGK primaries that have [Fe/H] 
measurements from SPOCS (IValenti &: Fischeii l2005l ). which adopts an uncertainty in [Fe/H] of 
±0.030 dex. The p roperties of these syste ms, which include two new proper motion-selected bi- 



nary systems (from iLepine &: Gaidoa l201lh , are listed in Table [TJ We are currently analyzing the 



remainder of our M dwarf sample, and will publish these results in a separate paper. 

We operated SpeX in the short cross-dispersed (SXD) mode with the 0.3x15" slit (2 pixels 
per resolution element), which yielded -R~2000 spectra spanning 0.8-2.4 /im. SpeX uses a standard 
nodding technique for sky subtraction. For our standard targets, we reached a minimum S/N~150 
per pixel, while for the calibrators, we obtained S/N~200 per pixel or greater. In clear conditions, 
we were able to achieve S/N~150 on an H=9 target with 10 minutes of integration. 



We extracted our spectra using the SpeXTool package (jCushing et al.ll2004l ). which automati- 
cally subtracts the AB nods, divides out flat-field exposures, wavelength calibrates, and performs 
optimal extraction. To achieve the best data quality, and avoid known instrument shifts with large 
telescope slews, we obtained flat-field exposures and argon lamp calibration sequences throughout 
each night. Each spectrum was telluric-corrected and flux-calibrated using a nearby AOV star, or a 
close spectral type if there was no nearby AOV star. Our telluric standards were always within 0.1 
airmass of the targets, and usually much closer. They were observed within 30-60 minutes of the 
target s. Telluric correct ion and flux calibration were performed with the xtellcor telluric correction 
code (|Vacca et al.ll2002l ^. 



3. Analysis and Results 

Our metallicity calibration is based on features in the H and -fC-bands. While our i^-band 
calibration is similar to RAlO, we use a different flux windows and continuum estimation regions, as 
described below. We also measured the equivalent widths (EWs) of two i?-band features, which we 
label the K i (1.52//m) and the Ca i features (1.62^m) based on their dominant elemental opacities. 
The decomposition of the ff-band features is shown in Figure [H The K i feature is weak and 
contains a telluric H2O feature. The -ff-band Ca i feature contains multiple atomic lines and is 
composed of the sum of two individual features, hereafter denoted Ca(l) and Ca(2). In order to 
perform our H-haxid analysis, we required S/N~150 per pixel, compared to the S/N~100 for the 
stronger K-haxid features. 
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3.1. Equivalent Width Calculation 

The first step in measuring the EWs was to define a standard width and location for each 
feature. The width had to be sufficient to contain the feature at aU strengths but not so wide as 
to contain nearby features. Our fixed feature widths, W, were: 





-band 


= 71 A, 


VFca,K- 


-band 


= 74 A, 


Ca(l),H- 


-band 


= 25 A, 


Ca(2),H- 


-band 


= 33 A, 




-band 


= 24 A. 



As discussed in N3.21 the feature widths were chosen to produce the best agreement between our 
analytic and Monte Carlo error estimates, signifying that the EW measurements were robust to 
small changes in the spectrum. These windows were positioned on centers C: 



CNa,K- 


-band 


= 2.2074 /im 


Cca.K- 


-band 


= 2.2638 fim 


Cca(l),H- 


band 


= 1.6159 /im 


Cca(2),H- 


-band 


= 1.6203 /im 


Ck,h- 


-band 


= 1.5171 ;um 



Calculation of the EWs also requires an estimate of the pseudo-continuum for normalization. 
To obtain this, we selected regions in each order that bracketed our features and contained no 
strong features themselves. The continuum regions are listed in Figure [TJ These regions were fit 
by a fourth-order Legendre polynomial using the IDL svdfit routine, and the fits were visually 
confirmed. 

We then accounted for radial velocity (RV), opting to use the strong features provided by the 
entire i('-band order and not the comparatively noisy ff-band. As a temp late, we used the h igh- 
S/N spectrum of the M1.5V star HD36395 from the IRTF spectral library (|Ravner et al-lbood ). A 
fit to the cross-correlation peak yielded a velocity offset, which was applied to all the feature and 
continuum regions for each target. 



EW measurement for each feature then proceeded in the standard fashion: 



EW, 



1 



AXi, 



(1) 



where the feature spans pixels i, Aj is the wavelength at pixel i, / is the intensity, and Ic is the 
psuedo-continuum intenstiy. We included in this sum fractional "edge" pixels which fell inside the 
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feature windows, weighted by the proportion of the pixel contained in the window. For the H-band 
Ca I feature, we combined the two components: 



EW 



Ca,H-band — EWca(l),H-band + EWca(2), H-band- 



(2) 



As suggested by RAIO, we accounted for temperature effects by including in our model H2O 
indicies in bo th the K arid iJ- bands, as the strength of H2O absorption in these bands tracks 
spectral type (jCovev et alj 120101 ). These indicies are defined as: 

(J"(2.180 - 2.200))/(J'(2.270 - 2.290)) 



H2O-K: 



H2O-H: 



(J^(2.270 - 2.290))/(J'(2.360 - 2.380)) 
(J"(1.595 - 1.615))/(J'(1.680 - 1.700)) 



(7'(1.680 - 1.700))/ (J'(1.760 - 1.780)) 

where {J-{a — h)) is the mean flux in the range bou nded by a and 6 in u m. A calibration of H2O-K 
based on the stars in the IRTF Spectral Library ( Ravner et al.ll2009l ) ] 
estimates listed in Table [TJ 



2009!) provides the spectral type 



3.2. Equivalent Width Errors 

We used two independent methods to e stimate uncertaintie s in ou r EW measurements. The 



first was an analytic formalism described in ISembach &: Savagd (Il992l ). Briefly, this method ac- 



counts both for errors in intensity and errors in the continuum fit, using the measured uncertainty 
from the extraction and the covariance matrix from the least squares fit to the pseudo-continuum. 
This yielded average relative errors for our EWs of 1-3%. 

Since this method did not account for errors induced by our RV correction, we also computed 
Monte Carlo errors. We performed 100 trials for each target, jittering each pixel by a draw from a 
normal distribution defined by its extracted error (~0.3-0.5%), and repeated the entire EW analysis 
on each iteration. The spread in the recovered EW values included errors in RV shifts, intensity, 
and pseudo-continuum fitting. In Table [H we report the more conservative Monte Carlo uncertain- 
ties. We selected our feature definitions by minimizing (through trial-and-error) the discrepancies 
between our analytic and Monte Carlo errors. For all except two metallicity calibrators, the Monte 
Carlo and analytic error estimates were consistent, their uncertainties being dominated by intensity 
and continuum-fitting uncertainties rather than RV uncertainties. 

We note that some of our observations were performed in conditions with highly variable cloud 
cover and water column. We diagnosed this issue by reducing individual frames of a selection of 
targets which had high S /N in each frame. In poor conditions, we found that the EW of the weak H- 
band K i feature could vary between frames by as much as 20-30%(~0.2A). We suspec t this is due to 



stron g nearby H2O features (as shown in Figure [T]). Simulating strong H2O features (jClough et al 



2OO5I ) in varying water column and calculating the resultant EWs showed that a realistic variation 



of 1mm of water column yielded AEW~0.2A, enough to account for the observed fluctuations. 
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3.3. Correlations 



We performed two least-squares linear regressions on our calibration targets, with metallicity 
as the response: one with the -fC-band features and the H2O-K index, and one with the i?-band 
features and the H2O-H index. For the K-hand, we found the best fit equation 



[Fe/H]K-band 
(T([Fe/H]K-band) 



0.132(EWNa) + 0.083(EWca 
-0.403(H2O-K) - 0.616, 
0.12. 



For the //-band, we found: 



[Fe/H]H-band 



0-([Fe/H]H-band) 



0.340(EWca) + 0.407(EWk) 
+0.436(H2O-H) - 1.485, 
0.12, 



where a is the dispersion of the residuals, and an average of the K and iif-band metallicity estimates 
also yields a residual dispersion of (T([Fe/H]avg) = 0.12. Residuals to these fits are shown in FigureEJ 
The observed fluctuations of the weak K i feature suggest that the H-hand metallicity estimate is 
more susceptible to atmospheric variability. For example, a variation of O.2A in the H2O feature 
would yield A[Fe/H]Ri0.08 dex. 

For comparison with other empirical methods, we present the residual mean squares {RMSp) 
and the a djusted squared multiple cor relation coefficients i?^ for our models and others (from 
RAlO and ISchlaufman &: Laughlinll2010l ). Our H and K-hand models provide equivalent fits: both 



have RMSp = 0.02 and they have = 0.74 and = 0.7 3, respectively. Our results are an 
improvement over those presented bv lRoias-Avala et al.l (j201l|) {RMSp = 0.02, Rl = 0.67), RAlO 
(R MSp = 0.02 R^ = 0.63), JA09 {RMSp = 0.04, Rl = 0.059), SLIO {RMS = 0.02, Rl = 0.49), 
and iNeves et al.l (|201ll ) {RMSp = 0.03, Rl = 0.43). Each quantity is reported for the calibration 
set with which it was created. 



As shown in Figure [21 our calibration is well-constrained for approximately —0.25 < [Fe/H] < 
0.3. It is poorly constrained for late (>M5), low-metallicity M dwarfs. We note that increasing 
the order of our model (e.g. to quadratic) does improve our fit, but this is mostly due to the single 
point with [Fe/H]spocs = —0.69, for which our calibrations overestimate the metallicity compared 
to the primary star (as does RAlO). 
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Discussion 



Figure [3] shows our calibration applied to our entire M dwarf sample, and demonstrates the 
equivalence of our K and if-band metallicity estimates, as well as the consistency of our imple- 
mentation of the RAIO K-hand calibration. We suspect that some of the scatter between the 
K and H-hand estimates is produced by the atmospheric variations discussed in §3.21 Other 
sources of scatter may include variations in stellar surface gravity, magnetic fields, and activity. 
Although progress is being made in characterizing activity-sensitive NIR features in M dwarfs (e.g. 
Schmidt et al.ll2012l ). our features have not been examined in thi s context. Howev er, we did es- 
timate the gravity sensitivity of our features using the BT-Settl (lAllard et al.l 120101 ) models. We 
applied our EW algorithm to stellar models with T=3100K, solar metallicity, and log(g) = 4.0, 4.5, 
and 5.0. This showed that the Ca features were not strongly gravity-sensitive, but that an increase 
of 0.1 in log(ff) yielded AEWNa~0.2A and AEWr^^O.iA. These changes translate to A[Fe/H]ss0.03 
and 0.04 dex, respectively. 



4.1. Calibration 

The main disadvantage of the H-hand method is that the i^-band features are much weaker 
than those in the K-hand, and therefore higher S/N spectra and non- variable conditions are required 
in order to accurately measure the feature strengths. However, the iJ-band offers three important 
advantages: 

• The if-band calibration allows an independent metallicity estimate from the X-band cali- 
bration, in a spectral region less plagued by thermal noise. The information content of the 
i7-band may be useful to other groups which are pursuing similar programs to ours. 

• The ff-band relation is superior to previous models in terms of R^. That this relation performs 
as well as it does despite its sensitivity to atmospheric effects is quite promising. Moreover, 
higher resolution observations can mitigate much of the atmospheric effects with detailed 
telluric modeling. 

• The ff-band model locates features that are strongly correlated with metallicity. These 
features and our calibration can be used as a starting point for M dwarf metallicity studies 
with high-resolution surveys like APOGEE. A recent analysis by Oil demonstrates that high 
resolution spectroscopic modeling of early M dwarfs is feasible in the J-band, raising the 
possibility that the i/-band might also be amenable to such techniques. 
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4.2. Metallicity of Planet Hosts 

We have also applied our calibration to several planet hosts, listed in Table [2l Our measure- 
ment of GJ 317 represents the first empirically calibrated spectroscopic metallicity estimate of this 
object. Our estimates are [ Fe/H]K = 0.26 and [Fe/Hjn = 0.31 dex, supporting the findings of 



Anglada-Escude et al.l (120111 ) and Oil that this star is metal-rich. 



We compare our results to other spectroscopic studies in Table El Generally, our results are 
consistent within errors with those of both RAIO and Oil. For the planets we have in common 
with RAID, our ET-band estimate differences have mean = —0.11 and a = 0.08 dex, and the -ff-band 
yields differences with mean = —0.14 and a = 0.12 dex. For planets we have in common with Oil, 
our i^-band estimate differences have mean = 0.04 and a = 0.09 dex and the -ff-band differences 
have mean = —0.01 and a = 0.10 dex. The most discrepant points, the H-hand estimates for 
GJ 849 and GJ 876, have high S/N which allowed us to probe frame-by-frame variability. The 
individual frames yielded highly variable (0.1-0.2A) EW measurements of the i^-band K i feature 
in both objects. 



4.3. Planet-Metallicity Relation 

Since we have obtained metallicity estimates for five M dwarfs with giant planets (of seven 
known) and four M dwarf planet hosts (of 12 known) without known giant planets, we can statisti- 
cally asses whether M dwarfs with giant planets are metal-rich. W e elect to use the nonparametric 



Mann- Whitney- Wilcoxon rank-sum test (jMann &: Whitnevlll947l ) on the averages of our K and H 



band metallicity estimates. This test determines whether one set of observations has preferentially 
higher values than another. For the samples of non-giant planet hosts (n = 4, median = —0.03) 
and giant planet hosts (n = 5, median = 0.27), we find that the distributions differed significantly 
(P < 0.01). The P- value is the probability of obtaining a rank-sum at least as extreme as the 
one obtained under the null hypothesis that neither group has preferentially higher metallicities, 
and our alternative hypothesis was that the giant planet hosts were more metal-rich. Thus, for the 
set of M dwarf planet hosts for which we have consistent metallicity estimates, we find that giant 
planet hosts are indeed preferentially metal-rich compared to M dwarfs with less massive planets. 



5. Conclusion 

We obtained NIR spectra of 22 M dwarf companions to FGK stars with known metallicities. 
From measurements of the strengths of two i?-band features and an H2O index, we have derived a 
new empirical metallicity relation for M dwarfs in the ff-band, parallel to that presented by RAlO 
in the X-band. This relation requires higher S/N spectra, but has smaller residuals. We used 
these calibrations to estimate metallicities for a subset of M dwarf planet hosts, producing the first 
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empirically calibrated spectroscopic metallicity estimate for the Jupiter-host GJ 317 and finding it 
to be metal-rich. We also used a nonparametric test to show that M dwarfs with giant planets are 
preferentially metal-rich. Finally, we anticipate that the metallicity-sensitive H-hand features used 
in our calibration will be critical starting points for metallicity analysis of high resolution i7-band 
spectra of M dwarfs from large surveys like APOGEE. 

We thank Kevin Luhman for obtaining observations. The authors are visiting astronomers at 
the Infrared Telescope Facility, operated by the University of Hawaii under Cooperative Agreement 
no. NNX-08AE38A with the National Aeronautics and Space Administration, Science Mission Di- 
rectorate, Planetary Astronomy Program. We acknowledge support from the NAI and PS ARC, and 
NSF grants AST-1006676 and AST-1126413. The Center for Exoplanets and Habitable Worlds is 
supported by the Pennsylvania State University, the Eberly College of Science, and the Pennsylva- 
nia Space Grant Consortium. The authors acknowledge the significant cultural role and reverence 
that the summit of Mauna Kea has within the indigenous Hawaiian community. We thank the 
anonymous referee for valuable comments and suggestions. 
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Table 1. M Dwarf Metallicity Calibrators 
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-0.13 
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GJ 872 A 
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3.18±0.05 
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Note. — 



The list of M dwarf companions to higher mass primaries which are used as calibration targets for the metallicity relations. 



'^Spectral types are estimated from the H2O-K index value, following ICovev at alj (|2010r i. 
Primary metallicities from SPOCS, which adopts an uncertainty in [Fe/H] of 0.03 dex jValenti fc FischeilboOSi) . 



"The 2MASS IDs are provided for these targets, which are new proper motion-selected binaries from lLepine fc GaidosI (|2011 ). 
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Table 2. M Dwarf Planet Hosts 



Name 


[Fe/H]K 


[Fe/H]H 


RAIO 


Oil 


Known Giant Planet Host 


GJ 436 


-0.00 


+0.02 


-0.00 


+0.08 ± 0.05 


no 


GJ 581 


-0.04 


-0.09 


-0.02 


-0.15 + 0.03 


no 


HIP 79431 


+0.53 


+0.50 


+0.60 




yes 


GJ 649 


-0.07 


-0.05 


+0.14 




no 


HIP 57050 


-0.03 


+0.05 


+0.12 




no 


GJ 849 


+0.31 


+0.22 


+0.49 


+0.35 + 0.10 


yes 


GJ 876 


+0.26 


+0.11 


+0.43 


+0.12 + 0.15 


yes 


GJ 179 


+0.20 


+0.18 






yes 


GJ 317 


+0.26 


+0.31 




+0.20 + 0.14 


yes 



Note. — Metallicity estimates for the observed planet hosts, along with estimates from 
Oil and RAIO. Errors for both the H and K-band metallicity estimates are a = 0.12 dex. 
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Fig. 1. — (A,B) The H and X-band orders of GJ 436, with continuum regions marked in hght gray, 
feature regions marked in dark gray, and the continuum fit overlaid. The left and right limits of the 
unshifted continuum regions are noted in /im, and the spectra have been normalized by a constant. 
(C,D) High-resolution, continuum-normalized spectra of the //-band Ca i and K i features in G J 



436. Overlaid are: a model spectrum produced using l inelists from VALD (IPiskunov et al 



updated log{gf) values from lMelendez &: Barbuvl (jl999l ). and the SYNTH3 c ode (iKochukhov 



1995 ) 



2003) 



and convolved to i?~30,000; an observe d spectrum of GJ 4 36 at iZ~30,000 (jBender et al.ll2005l ): a 
high resolution telluric model spectrum ( Clough et al. 2005 1) that is dominated by narrow CO lines 
in the Ca i feature and by H2O in the K i feature; and our IRTF observations of GJ 436 (i?~2000) 
with the continuum shown and offset by 0.2. The white regions show the feature windows. 



- 16 - 



0.4 
0.2 

0.0 

-0.2 
-0.4 



I ^ 

D 



6- 

5- 

J 4- 

3 - 
2 



□ 
o 



« ° DO 

On o 



9 



□ 



□ □ 



K-Band □ 
H-Band ^ 



□ 
❖ 



M7 
M5 
M3 
MO 



-0.8 -0.6 -0.4 -0.2 0.0 0.2 0.4 

[Fe/H] 



0.6 



Fig. 2. — Top: The r esiduals between the K/H band metalhcity estimates and the SPOCS 
(jValenti &: Fisched l2005l ) primary metalhcity. Bottom: The distribution of calibrators in color 
and metalhcity, using visual and 2MASS J magnitude (jSkrutskie et al.ll2006l ). showing the well- 
constrained regions for the calibrations: approximately —0.25 < [Fe/H] < 0.3 and MO to M5. 
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[Fe/H]K 

Fig. 3. — Top: Our A'-band mctallicity estimates compared to the RAIO estimates for the metal- 
Ucity caUbrators we have in common. In both panels, the hne represents equivalence between the 
two estimates. The results are consistent. Bottom: A comparison of the K and H-hand metallicity 
estimates for all program stars. The i^-band estimates track the X-band estimates well. GJ 317, 
GJ 849, and GJ 876 are marked. 



